Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



ENERGY INJECTION IN GRB AFTERGLOW MODELS 

GUDLAUGUR JOHANNESSON, GUNNLAUGUR B.IORNSSON AND ElNAR H. GUDMUNDSSON* 

Draft version February 5, 2008 

ABSTRACT 

We extend the standard fireball model, widely used to interpret gamma-ray burst (GRB) after- 
glow light curves, to include energy injections, and apply the model to the afterglow light curves of 
GRB 990510, GRB 000301C and GRB 010222. We show that discrete energy injections can cause 
temporal variations in the optical light curves and present fits to the light curves of GRB 000301C as 
an example. A continuous injection may be required to interpret other bursts such as GRB 010222. 
The extended model accounts reasonably well for the observations in all bands ranging from A-rays 
. to radio wavelengths. In some cases, the radio light curves indicate that additional model ingredients 

. ' may be needed. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanisms: non-thermal 
— shock waves 
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1. INTRODUCTION 

After a quarter of a century of modest progress, a breakthrough in gamma-ra y burst (GRB) research came in 1997, 
when X-ray and optical afterglow emission was first detected from GRB 970228 l|van Paradiis et aT1ll997HCosta et alJ 
' 119971) . showing a redshift of z = 0.695. This confirmed the cosmological origin of GRBs and to date, cosmological 
Q\ ] redshifts have been obtained for about 50 afterglows 2 . 

£N) ■ Initially, afte rglow light c urves were often sparsely sampled and had a natural explanation in the standard fireball 
' model (see e.g. iPiranl 120(151 for a recent review). Many showed a smooth power-law decay followed by a steepening 
one to two days after the burst. This tempor al behavior has been interpreted as being due to a relativistically 
^sO . expa nding collimated outflow fe.g. lRhoadslll999l hereafter R99) or more recently to a structured jet viewed off axis 
O ■ (e.g. lZhang fc Meszarosll2002t iRossi et al.ll2002fl 

The standard fireball model assumes a single energy release event with subsequent interaction with an ambient 
O ■ medium with constant density or a stellar wind density profile. It was also generally assumed that the electron energy 
distribution index, p, was independent of tim e and had values between 2-3, based on theoretical arguments and energy 
considerations (e.g. iZhang fc Meszarosll2004l) . 
+2 • It is quite remarkable that the light curves of most observed afterglows where sufficiently detai led data coverage 
' exists, can be e xplained within the standard fireball model or slight improvements thereof (e.g. iPanaitescul feOOS: 
. . \ iZeh et ail] 12006). This most likely reflects the underlying universality of the afterglow generation mechanism. For 
^ . some bursts, however, interpretations with the unmodified standard model, have proved difficult due to bumps in the 
light curves or because some basic model parameters turn out to have values outside the range assumed in the model. 
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In particular, the slo pe of the power law energy d istribution of t he ra diati ng electrons, p, is fou nd to be less than 
2 in some cases (e.g. iPanaitesc ni fc Kumarll2001aj) . iDai fc Chengl lj2001fl and B hattacharval ((20011 have explored the 
light curves and spectral signatures in models with y < 2. Numerica l calculations of relativistic shock acceleration 
suggest a uni versal value of p sa 2.2 — 2 .3 (e.e. lAchterberg et alJl200lj) . Other work s uggests an even higher value of 
p « 2.6 — 2.7 l)Lemoine fc Revenull2006l) . This has been referred to as the p-problem llBj^ms^on etaL| |2002() . Recent 
afterglow observations do indicate, however, that the value of p may not be universal l|Sehn et al.ll2005(K 

The p-problem arises when the standard model is being used to interpret light curves that decay much slower than 
sugge sted by the model. Such a slow decay had in fact been anticipated ( Rees & Mes zaroslll998t iSari fc Meszarosl 
2000, hereafter RM98 and SM00, respectively) by the natural assumption that the mass ejected in each event (and 
the associated energy) is released with a range of Lorentz-factors rather than a single value. The distributed energy 
injection into the fireball sustains higher fluxes for longer periods of time, thus resulting in slowly decaying light curves, 
without requiring p < 2. The cause of the slow decay is thus shifted from the electron population to the properties 
of t he source itself. This extension of the standard m odel was, however, not used to interpret actual data until later 
(e.g. iNakar fc Piranll2003l iBiornsson et alJl2002L [2004) . As continuous energy injection can sustain slower light curve 
decay rates with p > 2, it is important to explore fully the properties and light curve signatures of such models. 

I n addition to the above q uoted references, energy injection in the context of GRBs has previously been considered by 
e.g. lPanaitescu et al.l l|1998() who give an example of how the afterg low of GRB 970508 ma y be explained by refreshed 
shocks, although a spherically symmetric model may also wor k. iC ohen fc Piranl II1999T) explore the self-similarity 
of spherical blast waves with continuous energy injection, while IPailTTuri)2000l 120011 ) c onsider a model in which a 
pulsar is assumed to provide the additional energy, a variant of which was also considered by Zhang & Meszaros ( 2001) . 

* Science Institute, University of Iceland, Dunhaga 3, IS— 107 Reykjavik, Iceland, e-mail:gudlaugu, gulli, einar@raunvis.hi.is 
2 See http://www.mpe.mpg.de/~jcg/grbgen.html for the most recent list. 
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iKumar fc PirarJ (J2000) explore a toy model of the in teraction of two colliding shells and estimate the flux contribution 
from both the forward and the reverse shocks. Then. IZhang fc Meszarosl 1)20021) consider the intera ction of two colliding 
shells in detail and determine a criteria for such collisions to be mild or vi olent. |Panaite scul ll2005|) considers additiona l 
energy injection in interpreting observations of several GRBs. In addition. iGranot et al.l<)2003[) and lHuang et alJ l)2006f) 
conclude that discrete energy injections are needed to properly account for the optical afterglow of GRB 030329. All 
adopt different degrees of approximations on various model ingredients. 

In recent months, observations by the Swift satellite have shown that the early light curve behavior may not be as 
smooth as previously thought. The origin of these variation s may be due to the tra nsition from the prompt phase to 
the afterglow phase or due to a strong initial reverse shock llPanaitescu et alJ| 2006D. The slower decay ing portions of 
the light curv es may be due to contin uous energy injection l)Nousek et al.ll2005HGranot fc Kumar 2006). In particular , 
GRB 050319 (|C usumano et alJl2006[) and GRB 050315, have shown this kind of a behavior fsee Tcihncarini et all l2005. 
for the latter burst and further examples). 

In this paper we present a detailed analysis of relativistic fireballs with either discrete or continuous energy injection. 
We show that this has interesting observational consequences in addition to the slower d ecay and provid e s clea r 
signatures in the light curves. An example of these effects has already been presented bv iBiornsson et aP ([2004), 
where a discussion of the polarization properties can also be found. The clearest of those signatures are bumps in 
the light curves due to refreshed shocks, which may in fact have been observed in e.g. GRB 021004 and GRB 030329. 
These signatures provide additional information on the hydrodynamic properties of the fireball. We give examples of 
fits to three GRB afterglows and show that even with the additional energy injection, there remains some discrepancy 
between the model and the data, especially at radio wavelengths. In our approach we try to keep the treatment as 
self-consistent as possible and improve on several of the approximations of some of the previous work. 

The paper is structured as follows: In Section |3 we review the effects that energy injections have on the dynamics 
of the fireball and outline our approach to the problem. In section |3 we present approximations for simple modeling 
of the dynamical evolution. In Section 0] we discuss the ingredients in the radiative component of our model and in 
Section[5]we give examples of bursts that may be interpreted with this modified model. Finally, in Section|H|we discuss 
our results and conclude the paper. 

2. THE DYNAMICS 

We begin this section by considering energy and momentum conservation. We follow the approach of RM98 and 
SM00 and assume that energy is released over a period that is short compared to the duration of the afterglow. 
We then specialize the method to three cases : An instantaneous energy injection in Section 12. II for easy comparison 
with R99 (see also iPaczvnski fc Rhoadslll993[l . several discrete injection events in Section l2~2l and finally continuous 
injection with a power law distribution of the Lorentz factors in Section 12.31 We assume the dynamical evolution to 
be adiabatic in all cases. 

Assume that the total ejected mass with initial Lorentz factor higher than T can be written as M(> T) = Ma + AM, 
where Mq is the mass ejected with the highest Lorentz factor To, and AM ~ j T ^(dM / dT)dT , where dM/dT is a 
function specifying the distribution of ejected mass with Lorentz factor T. We prefer to use x = T/T as a variable, 
and introduce a dimensionless function F(x), defined by F(x) = M(> T)/Mo, in the interval 1 > x > r m /ro, where 
r m is the minimum Lorentz factor occurring in the ejected mass. We then have dF/dx — (To/Mo)dM/dT. Note that 
in general dF/dx < 0, since in realistic scenarios we must have AM > 0. Note also that F(l) — 1. The original fireball 
scenario with an instantaneous energy release is equivalent to dF/dx = 0. 

The energy of a mass element dM is dE = TdMc 2 , where c is the velocity of light, and the total energy of the ejecta 
with Lorentz factor greater than T is therefore given by 

E(> T)=E + AE = E + j dT = E Q (l+ j x^-dx\ , (1) 



where Eq = TqMqc . Integrating by parts we find that 



E = E ]xF(x) + j F(x)dxj . (2) 

The expanding shell decelerates by sweeping up the ambient medium. Assuming that the instantaneous total mass of 
the swept up medium is M e , the total energy of the shell can be expressed as E = T(-fi[M(> F) + M e ]c 2 ), where ji 
is the average internal Lorentz factor of the particles in the expanding shell. 3 Using this and equation (0) above, we 
find that the conservation of energy may be written as 

r 7i (F(z) + /) = T [xF(x) + h] + /, (3) 

with / = M e /M and 

h = j F(x)dx. (4) 



3 To be fully consistent, one would also need to distinguish between the internal Lorentz factors of the ejected shells and that of the 
shock-front. This distinction does not, however, change the final results as the total energy of the shell cancels out in the derivation of 
equation JSJ below. 
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Note that the swept up mass ratio, /, increases as T decreases. 

In the preceding equations we have assumed that the slower moving mass shells ejected by the central engine catch 
up with the decelerating shock front as soon as the Lorentz factors of the two coincide. This is a simplification, as the 
speed difference between the ejected mass elements will cause a time delay. The dynamical evolution is however not 
affected by this, but it should be kept in mind when considering reverse shocks originating from shells colliding with 
the shock front. A method for including this time delay will be discussed in section l2~2l 

Before impact with the expanding shock front, each mass element has momentum dP = TdMftc, where (3c is the 
velocity of the mass element in the burster frame. Hence the total momentum before impact is 



P = P + / dP = P - M c / {T 2 x 2 - 1) 1/2 — rfx, 
Jr Jx ax 



(5) 

with P = T Mq(3oC. As the momentum of the hot shell is P = T^i\M(> Y) + M e ](3c, momentum conservation results 
in 

T li (3{F{x) + f)=T (3 -I 2 , (6) 



where 

/"'■ ,IF 

(7) 



Dividing equation @ by equation ©, we obtain an expression for the dimensionless velocity of the expanding fireball, 

a = r oA> - J 2 ,„x 

' T [xF(x)+h]+f [) 
We then find that the Lorentz factor is given by 

r i r [ X F(x) +h] + f 



Vl^W V[To[xF(x)+I 1 }+f}^[r f3o-I: 
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This is an implicit equation for T, and describes the evolution of the Lorentz factor with radius (or time), once the 
distribution of the injected Lorentz factors, F(x), has been specified along with the initial conditions of the fireball 
(e.g. To and Mo). A description of the swept up mass, /, is also required (see eq. [TU] below). Equation is the key 
equation for obtaining the dynamical evolution of the fireball and has to be solved numerically in all realistic cases. 
Since we take the evolution to be adiabatic, the dynamics can be solved for separately and the radiative properties 
computed afterwards. 

For a complete description of the fireball evolution we must solve four differential equations (R99) 

dr df _ 1 dt _ l + z df _ ftp 2 

dh^^ ' dT b ~f' df b ~ L 2 (l + /3)' ^~Mo r ' ( } 

Here, r is the radial coordinate in the burster frame and t b , t' and t are the burster frame time, the comoving frame 
time and the observer frame time, respectively. The solid angle into which the energy is beamed is estimated by 
SI = 27r(l — cos(# + c s t'/(ctb))) where #o is the initial opening angle of the collimated outflow and c s w c/y/3 is the 
sound speed in the comoving frame, assumed to be relativistic. The mass density of the external medium is p, which 
in general is a function of radius. 

2.1. Instantaneous Energy Injection 

In the standard fireball scenario it is assumed that the energy is released instantaneously with x = 1 and F(x) = 1. 
We then find that xF(x) + Ii = 1; 1% = and equation J^J reduces to 

r + / r + / 



v/rr + /] 2 - [r§ - l] v/i + 2r / + /2' 
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which agrees with the result of lPaczvnski &: Rhoadsl l|1993D and is the starting point of the R99 analy s is (see his eq. [4]). 
Some authors have based their analysis of GRB afterglows on this equation (e.g. lSalmons'onll2003t IBianco fc Ru ffini 
12001 . 

GRB 990510 is an example of a burst where the afterglow may be readily interpreted with a single instantaneous 
energy release event. We present our fit to the observations in Section l5~Tl 

2.2. Discrete Energy Injection 

Assume that in addition to the initial injection Mo, the central engine expels several shells with masses Mj and 
Lorentz factors Tj. Converting to dimensionless variables, Xi = T^/To and m, = M^/Mq we have that 

F{x) = l + ^m i H(x i -x), (12) 
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where 



is the Heaviside step function. We then obtain, xF(x)+Ii = XiiriiH(xi—x) and I? = — X] m iV^o x i lH(xi—x). 

Note that by the assumptions of our approach discussed in the first two paragraphs of Section [3 the Xi 's form a 
decreasing sequence. Inserting equation ^| into equation 10 wc obtain 

r, _ r (1 + J2 x l ra i H{x l - x)) + / 



[T (1 + J2x imi H{ Xl - x)) + ix - [r /3 + E^v /r §^- 1 ^(^ - ; 

Because the Heaviside function is zero until Xi > x, it is clear from this equation that the dynamical evolution between 
injections is similar to the instantaneous case as the evolution of T depends only on / in the same way. 

When a shell collides with the sh ock front, a reverse shock may propagate back through the shell if the collision 
is violent l|Zhang fc Meszaros1l2002|) . Depending on the circumstance s, radiation fro m the reverse shock may make 
a substantial contribution to the afterglow light l)Nakar fc Piranir2004|) . According to lZhang fc Meszarosl l)2002f) . the 
liklihood of getting a reverse shock depends on the relative energy and speed between the shells and the width of the 
incoming shell. As the energy is known and we neglect the dynamics in the shell collision, only the relative speed is 
needed. Assuming that only the shock front slows down due to interaction with the surrounding medium and that 
the trailing shell travels with constant speed, shells ejected with the average Lorentz factor of the shock, r a , will have 
reached the shock front at r in time t. R e placin g x with x a = r o /To in equation @ is then sufficient to account for 
the delay. As shown bv iKumar fc Piranl l)2000|) . the Lorentz factor of the inco ming shell is r n ss 2T in a constant 
density environment for most of the evolution. Combined with the results from Zhang & Meszaros ( 2002) (see their 
figure 3) this indicates that for a reverse shock to result, the energy of the incoming shell must be at least twice the 
energy of the shock front. 

We emphasize that in our numerical calculations we assume the shell collisions to be instantaneous and we neglect 
the effects of the interaction on the light cu rves. These effects are expected to be small considering that the relative 
Lorentz factor of the collision is T r m 1.25 (jKumar fc Piranll200fl) . We estimate that the shells would have to be at 
least several hundred light seconds thick for the effects of the interaction to be measurable with current instruments; 
these would also be strongly smoothed because of the equal arrival time surface (EATS). 

An example of an afterglow that may be interpreted with discrete energy injections is that of GRB 021004. 
iBicirnsson et alJ l|2004D showed that its properties, including polarization, could be explained by using such a model. 
Energy injection does affect the polarization levels calculated using the model of Ghisclli ni fc Lazzatil l)1999j) as these 
are most strongly determined by the Lorentz factor. Changes in the Lorentz factor therefore manifests itself not only 
in the light c urve but also i n pol arization data. An analysis of a more detailed data sets for GRB 021004 has been 
presented by Postigo et alJ l)2005|) . The late time afterglow and the host galaxy of this burst have been studied by 
iFvnbo et all J2005D . 

The afterglow of GRB 030329 also showed variations b oth in the optical a nd the polarization light curves. These 
have previousl y been interpreted with energy inject ions ijGranot et al.ll2003[) . although other alternatives have also 
been explored ijBerger et al.ll2003l iHuang et al.ll2006() . An application of our approach to that particular burst will be 
presented elsewhere (Guizy et al. in preparation). 

2.3. Continuous Injection 

Here, we follow the approach of RM98 and SM00 and assume that F(x) = x~ s , with ,s > 1, i.e. the energy is 
dominated by the lowest Lorentz factors. (The case s < 1 is equivalent to instantaneous injection as the energy is then 
dominated by the highest Lorentz factors). Then xF(x) + 1± = Z/Tq + 1, where 

Z=^- i (x-^-l). (15) 



To evaluate I2, we assume that ^/TqX 2 — 1 « Tqx, which is always a good approximation if 3> 1. We then find 
1% ~ —Z. The dynamical equation for T (eq. |S]) in this case becomes 

r + f + z 



y/l + 2T Q Z(1 - fa) + 2r /(l + Z/To) + P 

Note that for Z = 0, we recover the result for the instantaneous case (eq. ^3)- Also note that the second term 
under the square root is ~ when To S> 1. When considering models with continuous injection, we numerically solve 
equation (|16fl . since in all relevant cases 3> 1. 

3. SIMPLE APPROXIMATIONS 

The quantities that determine the dynamical evolution of the fireball are the dimensionless mass ratios F(x) and /, 
representing the injected energy and the swept up mass, respectively (see eq. F(x) is an input quantity in the 
model, and hence each case has to be solved separately. On the other hand, / is only dependent on the evolution of 
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r, p and Q. In general, p oc 7'~ 9 , where g = for constant density interstellar medium (ISM) and g — 2 for a wind like 
environment. The evolution of / can easily be obtained in two cases: the confined outflow case when c s t' <C 9qv and Q 
may be considered to be a constant. In this case the dynamics is that of a spherical outflow. The second case is that of 
a laterally spreading outflow, when c s t' ^> Q^r and f2 is large enough to neglect the initial opening angle. These cases 
were referred to as the power law regime and the exponential regime, respectively, in R99. That naming convention is 
however too limiting as we show below, since for a continuous energy injection there is no exponential regime. 
As a prescription for / in the confined case, we evaluate it as 

(3 - g)M 



while in the laterally spreading case we adopt (see R99) 

^-p{c s t'f. (18) 



dr M 

To obtain a full solution of the problem and the light curve behavior, these need to be consistently included in the 
calculation. 

App roximate solutions applicable in the instantaneous case have been discussed in detail by e.g. R99 and Sal monsonl 
(2003). We emphasize that the approximations can give useful insight into the properties of the fireball but they 
should, however, only be expected to be valid during limited time intervals in the afterglow evolution. This is partly 
due to the validity of the underlying assumptions in those limited intervals, and partly due to the influence of the 
EATS on the evolution (see Section ll~l"|) . Another factor is the spectral energy distribution (SED), that is not really a 
set of broken power laws (see Appendix^), but has an intrinsic curvature, and so do the lightcurves. In the following 
we first give a brief discussion of the discrete case but then concentrate on the case of continuous energy injection. 

3.1. Discrete injections 

If the energy is released in several discrete events (as discussed in Section l2~2l . the evolution after a shell impact can 
be approximated with the evolution following a single release event with the appropriate initial parameters. Of these, 
only the energy is relevant, being equal to the total energy that has accumulated in the shock front after the impact. 
Integration over the EATS then smooths out the transition. 

As mentioned in Section l2~2l the relative speed between the impacting shell and the decelerating shock front is an 
important factor in evaluating the reverse shock emission arising from the collision. The speed of the incoming shell 
equals the average speed of the shock front at the time of collision, and to estimate it, we need to know how the 
radius of the fireball, r, evolves with time in the burster frame, Using Y ~ yTo/2/, appropriate for a single release 
(see R99), and equations l|17fl and Q1U[1. it is easy to show that in the confined case the speed of the incoming shell is 
Fa ~ v(4 — g)Y. Although this should only be valid early on in the evolution, numerical calculations indicate that it 
is in fact appropriate during most of the evolution. 

3.2. Continuous energy injection 

All results derived in this section assume that the shock speed, We first note that as / is the swept up mass 

in units of the initial ejected mass, it is small compared to the other terms in equation (|16f) for a substantial fraction 
of the initial fireball evolution, but will dominate at later times. 

Early in the evolution of the fireball, when x is close to 1, Z < Yq and we recover equation (4) in R99, and an 
identical initial fireball evolution. The most interesting case, however, is when Z/Yo > 1. From the definition of Z we 

see that this occurs when x = Y/Yq < (s/(2s — l)) 1 '^ - , i.e. when x < 2/3 for s = 2 and for x < y / 3/5, when s = 3. 
In this case equation JHJ simplifies to 

r H , r ° + / + z . (19) 

V1 + 2Z/ + /2 

When both Yq and / can be neglected compared to Z, we have that Y w \J Z/2f. Dividing by To and solving for 
x = r/Fo, we find 



1 a \ 



2r n s - i 



Note that in the limit, s — > 1, we recover the power law behavior of the instantaneous injection case, where x oc / 1//2 
(R99). 

Confined outflow: Here, we adopt equation (|17J) . with flo m tt9q, as the solid angle into which the beam propagates. 
Inserting into equation lj2U|) . we find that 

1 s (3-g)A/ o y /(6+1) r _3/(s+i) ocr _( 3 - fl )/( s +i)_ / 21 n 



2Fo s — 1 f^op 



G 



TABLE 1 

Exponents for F v oc t~ a v~P for a laterally expanding beam. 





u m < u < u c 


V > u c 


a 
P 


3(3-s + 4/3)/(s + 5) 
(p-l)/2 


2(1 - s + 6/3)/(s + 5) 
p/2 



As in R99, we also evaluate the relevant times. From dt'/dr = l/(aTo2;), we find by integration that the time in the 
comoving frame is 

t > K ( _±+L- \ r (s+4- g )/(s+l) (22) 

Similarly, from dt/dtb — (1 + z)/(2cTqX 2 ), the time in the observer frame is given by 

* ( 3 + 1 
t oc 

\s + 7-2g 

Inverting this result and inserting into equation l|21|l . we find that T oc i _ ( 3 ~3 , )/( s + 7_2 3 , ) ) in agreement with the SMOO 
result. Also note that if we set s = 1 and g = 0, we recover the results of R99. 

Laterally spreading outflow: The treatment is similar to that of R99. The g — 2 case cannot, however, be solved 
analytically due to the coupling of the variables. We therefore only consider the case g — in this section (constant 
density ISM). The general case has to be solved numerically. 

We start from equation l|18[l . divide by dt'/dr and integrate to obtain 

/(3+.)/(i+.) 0C ^( i /3_ Ci)j (24) 
3 

in agreement with equation (14) of R99 if we set s = 1. Here, c\ is a constant of integration to be determined by the 
initial conditions. It becomes negligible when c s t' 3> 0$r. 

The exponential behavior discussed by R99 is seen to follow from equation (|24() . If s = 1 (formally equivalent to the 
instantaneous case), then equation (|24|) shows that / oc t' so that dt'/dr oc t' , and the time in the comoving frame 
increases exponentially with r. Therefore / also increases exponentially with r. 

On the other hand, if s > 1, the behavior becomes a pure power-law, / oc £' 3 ( 1+s )/( 2 + s ) anc j f x r (2+«)/(s-i) w hi cn 
shows that the exponential regime is a special case. It follows that / oc r 3 ( 1+s " r oc i( s_1 )/( 5 + s ) and T oc i~ 3 /( 5+s ). 

The physical reason for this change in behavior is that the dynamical evolution is now determined by the power-law 
energy injection rather than the lateral spreading. We note that in general, if the incoming shell has the same opening 
angle as the initial 9q, the resulting decay will be a mixture of a power-law decay for the part of the outflow within 
9q, and an exponential for that part of the flow that is outside of 0$. 

By assumption, the dynamical evolution is independent of the radiative properties of the fireball. To proceed, 
a p rescription o f the radiative processes is required. We adopt the standard synchrotron emission m echanism (see 
e.g. iPira nl 120051 for a review), that is widely assumed to be the source of radiation in GRBs (e.g. lTavanilfl996t 
Meszaros & R eeslll99'7tlSari et al.lll998j) . The details of the radiative component of our modeling are given in Section^] 
and in Appendix lAl 

To explore the approximate light curve behavior, F v oc t~ a v~v , we need to evaluate the characteristic frequencies 
v m , v c and the flux F Vm at the former (See Appendix lAl for definitions of these quantities). 

For the confined case, the temporal behavior of the characteristic frequencies, the F Vm , and the light curve decay 
indices, a, agree with those previously derived for the forward shock by SMOO (see Table 1 in their paper). 



For the laterally spreading case we find on the other hand 

roc^ 3/(s+5) , (25) 
roci («-i)/(.+B) j (26) 

^ m ocr 12 /( s+5 >, (27) 

i/eocr 2 '"" 1 )/^ 6 ), (28) 

F m oct (3s " 9)/(s+5) . (29) 

The corresponding light curve exponents are given in Tabled (Note that if we let s = 1 we recover the corresponding 
results for the instantaneous case). From the indices in Tableland the results of SMOO we see that if v m < v < u C) 
the steepening in the light curve when the observer starts to see the entire outflow surface is Aa = 012 — oil — 



24(2 -I- 0)/{{s + 5)(s + 7)), with (3 — (p — l)/2. For parameters typically inferred in afterglows (p — 2.5 and f3 ~ 0.75) 
and taking s = 2 as an example, we see that Aa ~ 1. When v > v c the result is Aa = 24(l + (3)/((s + 5)(s + 7)), with 
(3 = p/2 in this case. For typical parameters we now find Aa ~ 0.85. As the steepening of the light curve is similar 
in the two cases, the magnitude of the break does not indicate in which spectral region the observing frequency is. 
Clearly, a more detailed fit to the observations is needed. Inferring p < 2 from a standard model fit may be taken as 
a strong indication that a continuous injection may be involved. 



r ,( S +7-2g)/( S +l) 



(23) 



7 



4. LIGHT CURVES AND SPECTRA 

As is customary in GRB fireball models an d discussed in the previous section, we assume the radiation to be of 
synchrotron origin. We adopt the approach of IWu et alJ l)2004|) . although their method of calculating the radiation 
flux by integrating the synchrotron emissivity over the electron energy distributi on is however too c omputer intensive 
for our purposes. I nstead, we use a brok en power law approximation as in lWiiers fc Galamal l)1999l) , smoothly joined 
together simil ar tolGranqt fc Saril |2002). The detailed expressions are summarized in Appendix IAI 

Similarly to IWu et al.l l|2004f) 7we calculate the synchrotron self absorption coefficient directly and therefore the self 
absorption frequency is not included in our smoothly joined power law approximation. We also include inverse Cqmpto n 
scattering, but neglect higher orders as these are suppressed by the Klein-Nishina effect (Panaitcs cu fc Kumarl feOOCO. 
The Compton Y parameter is calculated from 

Y = t T e [ l 2 K(l)dj, (30) 



3' 



where Af e is the normalized electron distribution and r e is the optical depth to electron scattering. The modified 
cooling Lorentz factor is found from the implicit equation 

Jcic = y^y, (31) 

where, 7 C , is the usual cooling Lorentz factor (see ea. I A4|l . The Comptonized emission component is then found from 
the synchrotron spectrum using the scalings, v 1 ^ = 27^ m , v]P = 2~fcjc v c,ic and -Pmax = T e P max ijSari fc Esinl200l(l . 

4.1. The emitting region 

In the standard afterglow model the fireball is described by the self-similar formulation of Blandfo rd fc McKed l)1976f) 
(hereafter BM). Analytical sol utions are however only available if th e ambient medium has a density profile of the 
form n cx r~ 9 . Some authors, (|Nakar et alJ 120031 : lLazzati et al.ll2002fl . have adopted these solutions in models with 
modest density variations. Using their approach, we are able to reproduce their results only if we assume that the 
shock thickness follows the BM solution, A ~ r/T 2 . Care must be exercised here, since in this approximation the 
number of radia ting particles may no t be es timated correctly. To clarify, let us consider a fireball in a constant density 
ISM. Using the iBlandford fc McKed l)1976j) solution for the density behind the shock front and integrating over the 
swept up volume we find that the number of particles contained behind the shock is 

N p = ^ + 0(T-% (32) 

where r is the the radius of the swept up volume, and n is the (constant) ISM particle density. This result is appropriate 
for a homogeneous medium but leads to over- or underestimates, depending on the sign of the density gradient, in the 
number of radiating particles when used for variable ISM densities. In such cases it leads to incorrect flux estimates. We 
also remind the reader that density variations do affect the evolution of T which has to be calculated self-consistently 
to obtain correct fluxes. 

Instead of using the BM solution directly, we make the simplifying assumption that the shock front is homogeneous 
in the comoving frame. Using particle conservation and the jump conditions we find that the shock thickness in the 
burster frame is 

M e 

= 2Trm p r 2 (l-cos6)Tn /1 ^ 

appropriate for all density profiles. Here, m p , is the proton mass and, n' = 4Fn, is the density in the comoving frame. 
Note that whe n (1 — cosff) = 2, a nd M e = 4irr 3 n/3, we recover the BM shock thickness. Although apparently a 
crude estimate, iGranot et alJ l)1999l) showed that the difference between the calculated flux levels in such a thin shell 
approximation and the full solution are small. 

In all our flux calculations, we follow the evolution of the EATS. The total flux a rriving at the obser ver at a given 
time, t, is obtained by integrating over the EATS for a given observer frequency, v, l|Granot et alJll99fltl : 

FM = ^dfi J r 2( i-/3cos^ r dr dcos0 > (34) 

where, P', is the radiation power density in the comoving frame (see also appendix lA"|l . v' is the frequency in the 
comoving frame and 

t rcosO . . 

tern = ": ; I j (35) 

1 + Z c 

is the time of emission in the burster frame. Also, di, is the luminosity distance which depends on the cosmology. In 
our numerical calculations we assume a standard cosmological model with fl\ — 0.7, ft m = 0.3 and Hq = 70 km s _1 
Mpc" 1 . 



4.2. Energy injection effects on radiation properties 

To understand the effects energy injections have on the radiation properties of the afterglow it is best to start by 
considering the effects of a single injection. Since we neglect the dynamics of shell collisions, the Lorentz factor of 
the leading shell increases instantaneously as soon as another shell catches up with it. The subsequent evolution will 
then quickly relax to a solution which is similar to the one obtained if the total energy of both shells was injected 
instantaneously. As the energy of the shock front increases instantaneously, so does its luminosity. However, the 
increase in flux received by the observer will not be instantaneous but gradual because the integration over the EATS 
smooths out the trans ition. The t imescale of this gradual brightening depends highly on both the Lorentz factor and 
the jet opening angle l)Piranll2005[) . Energy injection in wider and slower moving jets will cause the flux to change on 
longer timescales and thus creates smoother bumps in the light curves. Shell interactions considered in more detail 
than in this work, would also contribute to the smoothening. 

5. EXAMPLES OF DALA MODELING 

In this Section we apply our model to observations and give three examples of how well our numerical results fit 
afterglow data. We use reduced x 2 deviation to assess the goodness of the fit. As our model extends to the entire 
electromagnetic spectrum, an un-weighted \ 2 is n °t adequate for our purposes. The number of points in the optical 
light curves often exceeds that in other wave bands by a factor of ten or more which tends to decrease the quality of 
the fit at the radio and X-ray wavelengths. To correct for this, we weight the contribution from each band with the 
square root of the number of data points in that band. We then normalize the result so that if the contribution from 
each point to the un-weighted \ 2 is unity, the weighted x 2 equals the total number of data points. The function we 
use to check the goodness of the fit is then 

2 1 A dp ^ 1 ^/ log(j^)-log(^)) \ 2 

Xd -°- f (A dp - aw) Ef=i U vm U \ ^ a + m,/F^ ) > [ > 

where N^p is the total number of data points, N\ c is the number of light curves, Ni is the number of data points for 
light curve i, A par is the number of parameters, Fij is the flux in light curve z, measured at time tj with the error 
AFij, and Fi(tj) is the flux from our model light curve i at time tj. Dividing by the number of degrees of freedom 
(d.o.f.) in the first factor of the equation provides the reduced \ 2 which we denote by x 2 f ■ The second factor is the 
normalization, and l/^/Nl in the sum over light-curves (lc) is the weight. Note that the fit is done for log(F) rather 
than F, as this is equivalent to a liner fit in magnitudes. 

Fitting the model calculations to the measurements is nontrivial since the number of model parameters can be 
large and the goodness of fit function, x d f , has multiple local minima. This makes local gradient based search 
methods unsuitable for our purposes. Instead we employ a global search method known as the Evolution Strategy 
(R unarsson &: Yad l2000). The essence of the algorithm is as follows: Initially, a population of randomly generated 
A search points (parameter sets) are generated within the specified parameter bounds and the initial variance of the 
normally distributed search distribution is set. These A points are then evaluated using the x d f function and sorted. 
The best pL of those are then used as parent search points and each one of them generates \ j \i new points within their 
Gaussian search distribution. However, before doing so the search distribution itself, or rather its standard deviation, 
is varied using the log-normal distribution. In this manner the most favorable search distribution is adapted to the 
function landscape resulting in faster search. Initially the search distribution reaches the entire search space. The 
evolutionary process is repeated for a number of iterations until the search distribution's standard deviation is below 
a pre-specified limit. 

In our fits, we did not allow the ^-parameter (the ambient density profile) to vary, but assumed it to be fixed at 
either g — (homogeneous ISM) or at g = 2 (wind). In the following sections, we set g — 0, as this provides a better 
fit in all three cases considered here. In addition, the angle that the line of sight makes with the jet axis (the viewing 
angle) is assumed to be zero. This is because it mainly affects the polarization properties of the afterglow that are not 
considered here. 

The complex landscape of the goodness function, caused by the non linearity and complexity of the underlying 
model, makes it difficult to find a unique parameter set that is a well defined global minimum of the goodness function. 
Degeneracy of the parameters can cause several minima to be considered as candidates for the global minimum. The 
1 — a errors we quote in our parameter determinations reflect the statistical uncertainty of the model parameters as 
determined from the measurements and do not include estimates accounting for uncertainties in the model (see also 
lPanaitescul2001l for a similar discussion). The results of our fits for three GRB afterglows are given in Tabled Several 
comments on the results are in order. As the afterglow light curves are weakly dependent on Tq, it is not constrained 
by the data. We only quote lower limits on it. The electron energy slope, p, is constrained both by the spectral slopes 
and the afterglow slopes and is the most constrained model parameter. In the case of GRB 000301C we only quote an 
upper limit on es, as it is poorly constrained, mostly due to the lack of A-ray data. Finally, we note that the I — a 
errors quoted in the table are generally not symmetric around the best fit values. 

5.1. GRB 990510 

GRB 990510 has well resolved aftergl ow light curves in many wavelengths and shows a smooth achromatic break 
during the second day of the afterglow ijHarrison et alJll99 9: Kuul kers et alJl2000t iStanek et al.lll99"9t iHolland et alJ 
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TABLE 2 

Values of parameters obtained by fitting the three afterglows. The errors quoted are the 1-<t estimates. 
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Fig. 1. — The afterglow light curves of GRB 990510 fitted with a model with an instantaneous energy rel ease. The m ulti-band fit is 
reasonably good, resulting in n f = 2.6 with the model parameters given in table_2] The X-ray data is from lKuulkers et alj 120001) . the 
optical data from Harrison ct al. (1999); Stanck ct al. (1999) and the radio data from Harrison ct al. 11999). The optical data is corrected 
for Galactic extinction of E(B - V) = 0.2 ISchlegel et alj|199l) . 

12000^ . This afterglow is a good ex ample of smooth lightcurves and can be interpreted within the standard fireball 
model ijPanaitescu fc Kumarll2001b(l . It is therefore an ideal burst to test the accuracy of our model and to compare 
its results to those obtained by others. 

Figure [I] shows the afterglow light curves in radio, optical and X-rays. The fit resulted in \\ o f = ^-6, quite an 
acceptable fit given the simplifying assumptions for the underlying dynamics and radiation and the fact we are fitting 
the entire electromagnetic spectrum wit h a single model. The best fit parameters are given in table Most of the m 
are similar to the parameters found bv iPanaitescu fc Kumarl l)2001bD . within a factor of two or three. The order of 
magnitude difference in e e , the fixed fraction of the total internal energy assumed to be in the electron population, 
stems from our different definition of 7 m . Although not shown here, the slope of the light curves vary continuously 
from a = 0.6 at 0.1 day to about a — 2.2 at 30 days. A rapid steepening sets in at 1 day after the burst but there 
is no sharp transition from a pre- to a post-break slope. We note that for this burst, 9 s» 1/r at 0.7 days, just before 
the steepening of the light curve. 

GRB 990510 was the first burst where a polarized afterglow was detected ijCovino et al.11199'^) . However, as there 
was only one positive detection it does not constrain the model very much, although it further strengthens the case for 
a synchrotron origin of the afterglow emission. This polarization measurement is easily accounted for by an off-axis 
viewing angle. It was not included in the fit presented here. 

5.2. GRB 000301 C 

The afterglow of GRB 00030 1C is well time resolved at optical and radio wavelengths and was observed in the near 
infrared and at mm wavelengths although no A-ray observations have been published. A short lived brightening was 
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Fig. 2. — A sample of the available afterglow light curves of GRB 000301C fitted with a jet profile and an additional energy injection, 
Ei = 2.1Eq at ti = 1.4 days after the burst (solid curves). The fit has Xd o f = ^-8 using all available data and is a reasonably good fit. 
The energy injection allows for both the sharp break in the optical/nIR bands and the smooth decline of the radio light-curve. The bump 
around 4 days seen in the R band causes some of the contribution to Xd f, a l° n S with the radio data. It rises and decays too sharply 
to be fitted with an energy injection. The model parameters are given in Table |5j The dotted curves show the best fit model without 
additional energy injection, resulti ng in Xn.n.f. = H- Se e text f° r further discussion. Optical and near infrared data are from I Jensen et all 
120011) :lRTioads fc Fruchteil 120011); IBhargavife UowsiH 12000): IMasetti eTaD 12000) and IB agar et all 120001) . Radio and mm data are from 
IBerger et al.1 120001 h lSmith et aTT 120011) and lFrail et alj J2003t) . All optical data is correcte d for Galactic extinction of E(B - V) = 0.053 
( Schlc gel et alj|19 98) and intrinsic SMC like extinction of Ag =0.1 using the profile of Pci 1 1992). A host contribution of 40 Hz 20 /^Jy has 
been subtracted from the 8.46 GHz data. 

observ ed in the op tical ligh t curves at abou t 3.5 days followed by a steep achromatic break around 7 days after the 
burst l | Sa gar et all 120001 IBerger et alJl2000|) . It has been suggested that this bump may be caused by micro-lensing 
ijGarnavich et alJ l2000). Their % 2 value of the fit is reduced by half by introducing the lens effect. 

Here, we fit the afterglow light curves with our model containing one additional energy injection, E\ at t\ = 1.4 
days post-burst. This is shown in figure [21 where a host contribution of 40 ± 20 /xJy at 8.5 GH z has been subtracte d 
from the data. The data in the figure is corrected for Galactic exti nction of E (B — V) = 0.053 (Schle gei et alj fl9981. 
and an intrin sic SMC like extinct ion of Ab =0.1 using the profile of lPeil l|1992f) . This value for the intrinsic extinction 
was found bv lJensen et al.l lj2001|) which fitted a power law spectrum with various extinction curves to the measured 
SEDs. 

Fitting the data with a model with just an initial instantaneous energy release results in \d o f — H (dotted curves 
in fig. |2l ■ Adding the energy injection Ei, reduces this by a a factor of three to Xd o f = ^-8, and gives the parameters 
presented in table[3 Most of the reduction in the Xd o f ^ s due to the optical light curves, as the energy injection allows 
for a fairly constant emission at 8.5 GHz and the sharp steepening of the optical light curves at approximately 6-7 
days (solid curves in fig. El- The remaining contribution to the Xd o f ^ s due to the bump at 3.5 days in the optical light 
curves, that rises too sharply to be fitted with an injection event. This is because of the EATS and the homogeneity 
of the emitting region in our model, making the flux increase too smooth to fit the sharp rise and decline observed. 
This sharp bump may be due to an injection that is local in the sense that it only affects a small patch of the radiating 
shock front. We also note that the best fit is obtained by having the energy injection at ti — 1.4, i.e. just before 
the first optical data points, rather than at times just before the light curve bump. Such an early injection does fit 
the flux levels in the various bands better than an injection that is constrained to lie within the time window of the 
observations. 

Our model parameter values agree with those of IPanaitescul l|200H) within a factor of 2 or 3, except for quantities 
that are defined differently (such as e e ). Given the differences in approach and the different level of details in his work 
and ours, these differences are rather minor. 

Finally, although energy injection that is uniform across the forward shock front is unable to account for the bump in 
the lightcurves at 3.5 days, other explanations may apply besides micro-lensing. These could be density inhomogeneities 
in the ambient medium or an inhomogeneous injection that energizes the emitting surface only partly (patchy shell). 

5.3. GRB 010222 

GRB 010222 had well time resolved afterglow light curves at most wavelengths and can be interpreted with the 
standard broken power law in the optical bands. The shallow slopes of the light curves require a hard electron energy 
distribution (p < 2), if interpreted within the standard model. An unknown host extinction may influence the SED 
slope that in turn may lead to an inconsistent determination of p, between the spectra and the light-curves. 
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Fig. 3. — A sample of the available afterglow light curves of GRB 010222 fitted with a jet profile and a continuous energy injection (solid 
curves). The data has been corrected for Galactic extinction. The fit gave \\ f = 4-6, using all available data where the radio data gives 
most of the contri buti on. The mode l para meters are given in tablc|21 Optical data are from Galama ct al. (2003); Cowsik ct al. (2001); 
ISagar et aU 1200111 and lStanek et al] 1200111 . X-ray date from lin't Zand et alj 1200 J) and radio and millimeter data from lFrail et alj 120031) 
and references therein. The dotted curves show for comparison a fit with only instantaneous injection, resulting in Xd o f = 6- The increase 
in x\ o f ' s mos tly due to contributions from the radio light curve. 



iBidrnsson et alJ (|20Q2fl found that a continuous energy injection in an otherwise standard fireball model could 
explain both the light curves and SEDs without a hard electron energy distribution. We follow their suggestion and 
use continuous energy injection model to fit the afterglow data. The result is shown in figure |5] where the data points 
have been corrected for Galactic extinction of E(B — V) = 0.018 ilSchlegel et al.ll99S ) and intrinsic SMC like extinction 
of A B = 0.1 using the profil e oflPel lll99l. t hat was also adopted bv lGalama et all p00l|) . The host of GRB 010222 
is a bright starburst galaxy ijFrail et aHl2002]) and where available, the host contribution has been subtracted from the 
light curves. 

The fit resulted in Xd f = 4-6 with the model parameters presented in table [21 The fit is not particularly good, 
mainly because the model gives too high a flux at the radio wavelengths. Both the optical and the X-ray data are 
well accounted for. We note that the value of T TO is rather high (table O and the continuous energy injection ceases 
at approximately 0.1 day, which along with the jet geometry (r « 1/9 at about 3 days), accounts for the steepening 
in the light curves. If we fit the data without the continuous energy injection, we obtain Xd o f = ^-8 instead of 4.6, 
with the additional contribution mainly due to a worse fit at radio wavelengths. The model parameters, however, are 
similar in both cases; a narrow and energetic jet with a low magnetic field strength in a low density environment. The 
continuous energy injection thus better accounts for the radio light curves but otherwise the evolution in both cases is 
similar. We are therefore able to fit the optical data reasonably well within the standard model and still have p > 2. 

6. DISCUSSION 

We have presented a detailed description of a relativistic fireball afterglow model with discrete and continuous energy 
injection. We have given explicit examples of fits to observations for both types of injections and shown how they can 
provide a better fit to afterglow measurements than the standard model with single initial energy release. 

Although not a direct subject matter of the pre sent paper, we note tha t even polarization measurements can in some 
cases be accounted for with the extended model (Bio rnsson et al. 2004). However, interpreting afterglow polarimetry 
is not straightforward. First, polarization measurements are difficult as the polarization level in general seems to be 
small in afterglows. The data may therefore be sensitive to polarized contributions from other sources such as the host 
galaxy, which however is most likely to be constant. Second, our model is sensitive to the geometry of the outflow and 
the assumed homogeneity of the emission region. The expanding shock front is unlikely to be homogeneous during its 
entire evolution, especially with the added energy injections. This is because an injected shell will not expand laterally 
in exactly the same way as the original shock front. The transverse size of the incoming shell will then be smaller and 
its energy deposition in the shock front will be concentrated in the contact zone between the two. Accounting for the 
effects of the EATS, could then lead to a rise in the level of polarization shortly after the injection and even abrupt 
changes in the polarization angle if the center of the colliding shell does not coincide with the geometrical center of 
the front. 

The origin of the additional shells is most likely the same as the origin of the shells producing the gamma ray burst 
in internal shocks. It is highly unlikely that the matter ejected from the central engine always collides and collects into 
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one shell in the gamma ray burst event. The remaining shells will then be the ones to refresh the external shock later 
on in the evolution of the fireball. This might even be the case in most bursts, but goes unnoticed in many cases since 
the energy of the colliding shells must be of the same order or higher as the total accumulated energy of the front if 
their interaction is to produce a measurable brightening. 

When fitting GRB afterglow light curves with the standard fireball model, it is quite common to use power law 
approximations with a light curve break due to jet structure. The slope of the light curve, a, before and after the 
break is then used to determine the electron energy distribution index, p, which is then compared to the value of p 
inferred from the SED of the burst (see e.g. Sari et al. 1998). These results must, however, be carefully considered, 
as the value of a depends on s as well (see table [y, and may not be smoothly varying between two fixed asymptotes, 
although that seems to be the general case. This is especially true when the characteristic frequencies of the synchrotron 
spectrum, v m and v Cl are near the observing frequency, where the SED does not follow a power law. 

The light curve break time has traditionally been used to infer the jet opening angle from the assumption that it 
occurs when 9 ~ 1/r. This can lead to erroneous angle esti mates, because energy injections can delay the appearance 
of the light curve break. An example of this is GRB 021004 ijBiornsson et alJl2004|) . As a result, the opening angle will 
be overestimated and the energy co rrected for beaming may thus not be correct (iJo hanncss on et al.ll200ql . Studies 
based on such energy estimates fe.g. lFrail et al.ll200lHBloom et aLll2003l iGhirlanda et alJl2004j) . should include error 
analysis that accounts for this 'bias'. 
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APPENDIX 

SYNCHROTRON RADIATION MODEL 

We approximate the synchrotron radiation power density, P 1 , in the comoving frame (primed quantities) with a 
smoothly joined power-law at the characteristic frequencies. This takes the form 
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/ \ (1/3)(-ki) 
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for the slow cooling regime, v' m < v' c . Here p is the electron energy distribution index and v' m and v' c are the 
characteristic synchrotron frequencies corresponding, respectively, to the minimum Lorentz factor of the distribution, 
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and the Lorentz factor above which electrons are radiative, 
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(A3) 



(A4) 



Here, m p and m e are the mass of the proton and the electron, respectively, e e is the fraction of thermal energy 
in the electrons compared to the total thermal energy of the fireball, <jt is the Thompson scattering cross section. 
Furthermore, j3 is the dimensionless speed of the shock front and B' = J es87r(r — l)n' is the magnetic field strength, 
with €b the fraction of magnetic field energy compared to the total thermal energy of the fireball and n' = 4Tn is the 
density in the comoving frame. The characteristic frequencies are given by 

=*>!£• (A5) 

where i is either m or c. The maximum values for the power density are given by 

. , 2.234eVS' , , . 

PL.A = <f> P 5 (A6) 



11.17(p- 1) e-n'B 1 



(A7) 



max,s 3p_l m eC 2 

where, e, is the elementary charge and the parameters, <j> p and %p, are introduced here as in iWiiers fc Galamal (Jl999) 
to account for an isotropic distribution of angles between the electron velocity and the magnetic field. Comparing the 
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results from these approximations with the more accurate calculations using the approach from I Wu" eTall ll200l . we 
find the following p-dependence of the exponents and coefficients in the fast cooling regime to be 

K 1= 2.37-0.3p (A8) 

K 2 = 14.7 - 8.68p + 1. V (A9) 

P = 1.89-O.935p+O.17p 2 (A10) 

X P = 0.06 + 0.28p. (All) 

In the slow cooling regime, wc find similarly, 

K 3 = 6.94- 3.844p+0.62p 2 (A12) 

K 4 = 3.5-0.2p (A13) 

p = O.54 + O.O8p (A14) 

X P = 0.455 + 0.08p. (A15) 

The difference between these approximations and the more accurate calculations is generally less than 10% unless v c 
is close to v m . Then the difference can be up to 40%. However, for most parameter sets the characteristic frequencies 
are close, only very early in the afterglow evolution and this is therefore not a concern in the later evolution. Even 
then, the difference is mostly smoothed out because of the integration over the equal arrival time surface. We find the 
difference between the two solutions to be less than 10% for most realistic cases. We also remind the reader that the 
electron energy distribution is an approximation and little is known about the detailed microphysics of the shocked 
material. 
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